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I. INTRODUCTION 

A number of "two-state" models were proposed to 
rationalize the phase behavior of Poly(ethylene oxide) 
(PEO) and its solution thermodynamicsi^i^i'^i'^i^i^ Within 
these models the monomers are in dynamic equilibrium 
involving two interconverting states (Fig. The Flory- 
Huggins lattice and the mixing entropy of the chains are 
retained. Additional contributions are due to the mixing 
entropy of the different monomeric states and their inter- 
actions. While the two-state models were proposed for 
aqueous solutions of PEO they are also candidates for 
the description of other neutral water-soluble polymers 
such as Poly(N-vinylpyrrolidone) (PVP) and Poly(N- 
isopropylacrylamide) (PNIPAM)J: The equilibrium free 
energy obtained from these models can be expressed as^ 
F/kT = N-^(f)\-ii(j)+{l-(f)) \-ii{l-(f))+Xeff<f){l-<f)) where 
N is the polymerization degree. In distinction to the fa- 
miliar Flory free energy, the effective Flory interaction 
parameter Xef / is a function of both the monomer volume 
fraction, (j), and the temperature, T. In this formulation 
the specific features of a particular model are grouped 
into Xe//(T, (/)). In the following we consider the Self- 
Consistent Field (SCF) theory of brushes of "two-state 
polymers" in terms of XeffiT,(j)). A significant part of 
our discussion is devoted to brushes of polymers capable 
of undergoing a second type of phase separationi ^i^i^i^°i^^ 
Within a brush, this type of phase separation can lead to 
a vertical phase separation associated with a discontinu- 
ous concentration profileii2ii^ Our analysis focuses on the 
signatures of such phase separation. These include the 
non-monotonous variation of the brush thickness with the 
grafting density and the appearance of distinct regimes 
in the compression force profiles. 

This approach is of interest because of a number of 
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FIG. 1: Schematic picture of a brush of a "two-state poly- 
mer". Open and filled circles depict monomers in the differ- 
ent, interconverting, monomeric states. 



reasons: (i) Xeff{T,(j)) determines a number of impor- 
tant characteristics of the brush among them the con- 
centration profile, the distribution of free-ends and the 
force profile associated with the compression of the brush. 
Thus, a description of the brush behavior in terms of 
Xef / {T, (j)) accounts for the leading brush properties and 
facilitates the comparison of the predictions of the differ- 
ent models. The specific features of the individual models 
and their parameters come into play when the distribu- 
tion of the monomer states is of interest. However, as 
we shall discuss, even in this case it is convenient to first 
specify the brush characteristics in terms of Xef f (T, (p) . 
(ii) The formulation of the theory in terms of Xef / {T, 4>) 
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underlines the relationship to the measurable 

x{T,<p)^Xeff~il-<P)dxeff/d<f> 

as obtained from the study of the coUigative properties 
of the polymer solutionsii^ii^iiSiii x(T, 0) is helpful in 
determining the parameters of the models. In the con- 
text of brushes, the behavior of x{T, (f>) provides, as we 
shall discuss, a useful diagnostic for systems expected 
to exhibit a vertical phase separation within the brush, 
(iii) While our discussion focuses on the two-state mod- 
els cited above, the analysis can be extended to other 
modelsiSiiS that yield a (j) dependent Xeff- (iv) The SCF 
theory of brushes characterized by Xef f {T, 4>) suggests 
useful tests for the occurrence of vertical phase separa- 
tion. This is of interest, as we shall discuss in the final 
section, because of experimental indications that such be- 
havior occurs in brushes of PNIPAM. (v) The concentra- 
tion profiles obtained from the SCF theory are essentially 
identical to those derived^^ from the Pincus approxima- 
tion where the distribution of free-ends is assumed rather 
than derivedi22i2i In marked contrast, the compression 
force profiles are sensitive to the distribution of free-ends 
and the two methods yield different results. 

The two-state models differ in their identification of 
the interconverting states. Within the n-cluster modeiSiSi 
one is a bare monomer while the second is a monomer 
incorporated into a stable cluster of n monomers. In 
the remaining models one of the monomeric states is hy- 
drophilic and the other is hydrophobic. The hydrophilic 
state is preferred at low T while the hydrophobic state 
is favored at high T. In the Karlstrom modeli the 
two states differ in their dipole moment and their in- 
terconversion involves an internal rotation. The mod- 
els of Matsuyama and Tanaka,^ Bekiranov et ai^ and of 
Dormindotova** assume that the hydrophilic monomeric 
state forms a H-bond to a water molecule while the hy- 
drophobic state does not. The brush structure within 
the Karlstrom model was studied using numerical SCF 
theory of the Scheutjens-Fleer type^^'^'^ and allowed to 
rationalize the aggregation behavior of copolymers in- 
corporating PEG blocksJ^: The brush structure within 
the n-cluster model was studied using SCF theory^^ and 
by simulationsiS^ These reveal the possibility of a ver- 
tical phase separation within the brush giving rise to a 
discontinuity in the concentration profile. In turn, this 
was invoked in order to rationalize observations about 
the collapse of PNIPAM brushes The force profiles due 
to the compression of brushes described by the n-cluster 
model were also analyzed.^^2L28, ^j^g studies of brushes of 
"two-state polymers" focused on a particular model and 
were formulated in terms of the corresponding free en- 
ergy. This obscured common features between the dif- 
ferent models and hampered the comparison between 
them. For example, while a vertical phase separation 
is possible within all two-state models, this scenario was 
mainly studied for the n-cluster model thus creating a 
misleading impression about the physical origins of this 
phenomenon. 



Our discussion concerns a brush of flexible "two-state" 
chains, terminally grafted to a planar surface. We as- 
sume that the chains are monodisperse and that each 
chain incorporates N monomers. The surface area per 
chain, cr, is constant and the surface is assumed to be 
non-adsorbing for the two monomeric states. The free 
energy per lattice site is 

/„,(0,T)/fcr= (l-0)ln(l-0)-fXe//(<^,T)0(l-</.) (1) 

This form corresponds to the N oo limit. It is appro- 
priate for brushes of any N because the grafted chains 
lose their mobility and thus have no translational entropy. 
The application of our analysis to a particular case is il- 
lustrated for the Karlstrom model. However, most of our 
analysis is model independent in that Xef f [T, (j)) is not 
specified explicitly. The only assumption made is that 
x(r, 0) can be expanded in powers of 

where the Xi(T) are specific to a given model. For sim- 
plicity we further limit the discussion to systems where 
the first three terms provide an accurate description of 
x{T, (p). As we shall discuss, the presence of a third order 
term is the minimal condition for the possibility of a ver- 
tical phase separation within the brush. The power series 
expansion of x(^^j 4') is clearly related to the virial expan- 
sion of the osmotic pressure, tt. The two differ in that the 
second incorporates terms originating in the translational 
entropy. In discussions of the SCF theory foo{(t>)/kT was 
often approximated by the second and third terms in the 
virial expansion of the Flory-Huggins free energy. From 
this point of view it is important to note two points: (i) 
in the power series of x{T, (j)) all the coefficients are T 
dependent and can change sign, (ii) Use of x(r, 0) se- 
ries expansion including a (jy^ term corresponds to a virial 
expansion incorporating a cj)'^ term. 

The next five sections, HlllVIl are devoted to the 
model independent aspects of the SCF theory based 
on foo{(t>,T)/kT with dependent Xeff- In section Hll 
we formulate the analytical SCF model for Xeff{4>,T) 
as a generalization of the familiar SCF theory of 
brushes The concentration profiles and their 
moments are discussed in section Hill The technical de- 
tails corresponding to section IIIII are described in Ap- 
pendix ^ Section IIVI discusses the distribution of free 
ends while the technical details are given in Appendix 
IbI The force profiles associated with the compression of 
the brush are analyzed in section In every case we 
distinguish between brushes characterized by a contin- 
uos concentration profile and those exhibiting a vertical 
phase separation. Finally, in section IVll we illustrate the 
implementation of our analysis to the Karlstrom model. 
In particular, we obtain the corresponding Xef f{T, 4')j 
x(r, (j)) and the distribution of monomeric states in the 
brush. 
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II. THE SCF THEORY FOR XeffiT^cj)) 

Consider a brush of neutral and flexible polymers com- 
prising TV monomers of size a. Each chain is grafted by 
one end onto an impermeable, non-adsorbing, planar sur- 
face. The area per chain is denoted by a and H is the 
maximal height of the brush. Following refs. I3I3I the 
free energy per chain, Fchain, consists of two terms: an 
interaction free energy, Fint, and an clastic free energy, 
Feh Fchain = Fint+ F^i- The interaction free energy per 
chain is 



F 

J. nn 



kT 



H 



fooWdz 



(2) 



where the interaction free energy density fooi't') is given 
by I'^ strong stretching limit, when the chains are 
extended significantly with respect to Gaussian dimen- 
sions, the elastic free energy is^'^ 



kT 
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g{z')dz' I E{z,z')dz. 
Jo 



(3) 



Here E{z,z') = dz/dn characterizes the local chain 
stretching at height z when the free end is at height z' . 
g{z') specifies the height distribution of the free ends and 
obeys the normalization condition g{z')dz' = 1. 

The concentration of monomers, 4>{z), at height z is 
specified by 



a' f" g{z')dz' 
Since each chain consists of N monomers we have 



(4) 



H 



TV = / (j){z)dz 



(5) 



At the same time 



iV = 



dz 



E{z,z' 



(6) 



which can be regarded as a normalization condition for 
the function E(z,z'). 

The equilibrium 4>{z) in the brush is determined by the 
variation of the functional Fchain with respect to E{z^ z') 
and g{z') subject to the constraints (0) and jnj yielding 



E{z,z') = 



(7) 



and 



where B = 37r^/8A^^a^, A is the Lagrange multiplier asso- 
ciated with constraint and ^(0) — dfoo{(t>)/d(j) is the 
exchange chemical potential. Up to this point the SCF 
theory is identical to the familiar versions, as obtained 
for X = x{T). 

For dilute brushes, cr >> 1, immersed in a good solvent 
and when x is independent oi(j), Xeff{4>jT) = x(T) <C 
1/2, the chemical potential is linear in 4>, fJ-ifk) ^ . In 
this case, when binary interactions are dominant, eq. ^ 
leads to a parabolic concentration profilei^iiS At higher 
grafting densities, cr > 1, higher order terms become sig- 
nificant. These were typically handled by incorporation 
of the third virial termnSiiS However, as discussed in 
the introduction, deviations from these scenarios are ex- 
pected when Xeffi'f) varies with (j) and /i('/>), as obtained 
from (Q, assumes the form 



= -\n{l-4>)-l + Xeff{c^)-2xeSf 



(9) 



Since colligative measurements yield x(</') = Xeffi'P) ~ 
(1 — (j))dxeff {'(>)/ d(j:^'^ rather than Xefs{4') it useful to 
express as 



li{c^) - - ln(l -ct>)-l + XeffiO) - / Xiq^)dq> ~ q^xW 

(10) 

Finally, the Lagrange multiplier A is determined by the 
concentration at the outer edge of the brush, (jj^ = 4>{z = 
H) 



A = BH^ + fi{z^H) 

= BH^ - ln(l -cPh)-1 + XeffiO) - 

' " x{cP)dcP-(t>HX{<t>H), (11) 



In turn, of a free brush is set by the osmotic pressure 
at H that 
leading to 



at H that is, -Kosm{4>H) = 0^9 [/oo(0)/0] /d(l)\^^^^ = 



- ln(l - ct^H) -4^H- x{4>h)4>h - (12) 
In a good solvent (f)H — 0. 



III. THE CONCENTRATION PROFILES AND 
THEIR MOMENTS 

In order to obtain (f){z) it helpful to express (jHJ as 



(13) 



(8) 



where 




FIG. 2: The state diagram of a brush with = 
and A*' — 200 in the X2, l/o" plane. Vertical phe 
tion occurs in the hatched region. At • i^o = 
dashed line, at higher l/cr, corresponds to brushes 
an inflection point at altitudes that increase with 
boundary for lower 1/a corresponds to (f)- (f)+ coe: 
the grafting surface. 
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xi(t>)d(t> + (l>HXi(t>H 



4>H 



Equation Ijl^ll does not specify 
it yields 



Tip 



z) directly. 
Afi{(f))/B. The brush 



determined in terms of the monomer volume f: 
the surface, cj)Q — 4>{z — 0), leading to H = yO. 
and 



A^(0o) - A^(0) 
B 



4>{z) is determined by equation H15|l together 
normalization condition (|3J|, which relates 00 to 
ing density, l/cr. 

We now distinguish between two cases. In one the 
concentration profile is continuous while in the second a 
discontinuity occurs due to vertical phase separation. In 
the first case ^ may be expressed as 



dz 



(16) 



The concentration profile for a given a is fully specified 
by H15() and . A vertical phase separation in the brush 
results in a discontinuity at height Hf. At this altitude 
two phases coexist: a dense inner phase with a monomer 
volume fraction 4>^(Ht) and a dilute outer phase with 
(j)-{Ht). In this case the normalization condition (0) as- 
sumes the form 



>|^d0 (17) 
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FIG. 3: (j>{z) for different areas per chain a when x{4') ~ 1/2-1- 
X20^ (a) X2 = 1 (b) X2 = 1-05. In every case = 200.^*^ 



where 4>+{Ht) and (j)-{Ht) are determined by = 
IJ,{4>-) E^nd 7r((/)_|_) = 7r((/)„). 4>{z) is now determined by 
equation (|15|l together with the normalization condition 

03, 

It is of interest to consider the phase behavior when 
x((/>) is described by xl"/") = Xo + Xi0 + X20^- In the case 
of x{4>) = xo or x(0) = Xo + Xi't' the critical point, as 
specified by fo.{<t>)/d<t>'' = d^fooW/d^^ = 



(1 



d(t> ^ 



(18) 
(19) 



occurs at 0c = 0. This corresponds to the familiar case 
of a polymer rich phase in coexistence with a neat sol- 
vent. In this situation there is no vertical phase sepa- 
ration within the brush and the concentration profile is 
continuous. A second type of phase separation,.5*S*2iili 
associated with a discontinuous (j){z)^ is possible when 
higher order terms are involved. For x(0) = 1/2 + X20^ 
a critical point occurs at 0c = 1/2 and X2c = 1- In the 
vicinity of the critical point, for X2 1 and > 1/2 the 
coexistence curve is well approximated by the spinodal 
line aVcx>(0)/502 =0 
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FIG. 4: Plots of (j){z) (a) and g(z) (b) above and below the 
critical point for x(0) = 1/2 + X2<^^>^ cr = 17 and = 200. In 
all cases, the outer phase is swollen4- 



FIG. 5: Plots of (j){z) (a) and g{z) (b) for the case of two 
coexisting dense phases with x{4>) — 0.51 +X20^, u = 18 and 
N = 200i^ 



1 - 



1-4x2 







(20) 



leading to <j)± = ^ ± \\J\ — 1/X2- The state diagram 
of a brush in the X2 , 1/f plane when eq. H20() applies is 
shown in Fig. |21 Concentration profiles obtained from 
x((/)) of this form are depicted in Fig. El 

Experimentally, the brush thickness H is inaccessible. 
Certain experimental technique yield (\)(z) allowing one 
to obtain moments of (\)(z) 



/o^ z(\>{z)dz 



H 



z(t>(z)dz (21) 



jj^ z'^(j){z)dz 
j" (j){z)dz 



H 



z^(l){z)dz. (22) 



Other techniques, such as ellipsometry, measure {z)m^ 
As we shall discuss, the a dependence of the moments 
provides useful information on the brush structure. The 
details of the calculation of these moments are described 
in Appendix 1X1 

When (j^iz) is continuous, both moments increase 
smoothly with the grafting density. In marked contrast, 
vertical phase separation within the brush gives rise to 



a non-monotonic behavior. In particular, both (z) and 
(z^) exhibit a minimum at intermediate a. A vertical 
phase separation gives rise to a plateau in the H vs. a 
plot (Fig. IHl while in the plots of (z) and \/Jz^ vs. a 
it is associated with a minimum (Fig. |S1 and Fig. • 
The physical origin of this behavior is the partitioning 
of the monomers between the inner dense phase and the 
outer dilute one. The minima are traceable to the higher 
weight give to the inner phase. Since the inner phase 
is denser, the onset of vertical phase separation is asso- 
ciated with a decrease (z) and \/{z^- These features 
provide a useful diagnostic for the occurrence of a ver- 
tical phase separation in the brush. The SCF analysis 
in this section confirms earlier resultai^ obtained by uti- 
lizing the Pincus approximationiSSiSi As we shall discuss 
this is the case for properties that are insensitive to the 
precise form of g{z). In marked contrast, the compres- 
sion force profile (section V) does depend on g{z) and the 
SCF result differ from the one obtained from the Pincus 
approximation. 



IV. THE DISTRIBUTION OF FREE ENDS 

The SCF formalism allows to obtain the height dis- 
tribution of free ends, g{z). Current experimental tech- 
niques do not allow to probe g{z) directly. However, g(z) 
is of interest because it plays a role in the calculation of 
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FIG. 6: (z) as a function of a for different X2 values and 
N = 3004^ 



FIG. 7: ^ (z^) as a function of a for different X2 values and 
iV = 3004^ 



the compression force profile. When x = const the brush 
structure is dominated by the contributions of the second 
and third virial terms of Fmt- Three scenarios emerge. 
In a good solvent the ends are distributed throughout the 
brush and 17(2) is a smooth function vanishing at z = 
and z = H . When the brush is collapsed in a poor sol- 
vent the ends reside preferentially at the outer edge of 
the brush and g{z) diverges at H. In a 6* solvent g{z) in- 
creases smoothly with z but does not diverge i^Si^i As we 
shall see, a new scenario emerges when a vertical phase 
separation occurs. In particular, g{z) will then diverge 
at the phase boundary indicating localization of the ends 
at the boundary. We will obtain g{z) from the integral 
equation I^J . The details of the calculation are described 
in Appendix IbI 

g{z) of a brush with a continuous profile (Fig. ^) is 
specified by 



del)' 



y/Afi{cb) - A^((/.') 



(23) 



where 4> and z are related by . When a vertical sep- 
aration occurs within the brush, equation Q yields now 
two expressions, for 4>{z). At the outer edge, Ht < z < H 



a' f" g{z')dz' 
while at the inner dense phase, < z < Ht 



(24) 



(^) 



a 
a 



Ht 



g{z')dz' r" g{z')dz' 



E{z,z') a Jh, E{z,z') 



(25) 



In the outer region only free ends with Ht < z contribute 
while for the inner phase all free ends are involved. 



The expression for g{z) in the two regions are given 
below while the details of the derivation are presented in 
appendix IbI At the outer phase, Ht < z < H 



g{z) = z- 



while in the inner phase, < z < Ht 



(26) 



\Ht)~^-{Ht) , (l)H 



d(j)' 



VH^ - z^ 
B 



Hf - z2 X 



d<j)" 

0, T^Mpp^Mpy 



(27) 



The first integral allows for the contribution of the inner 
phase and the second for the contribution of the outer 
phase. g{z) (|27|) at the interval < z < Ht diverges at 
the phase boundary z — Ht- g{z) (|26() at the interval 
Ht < z < H diverges at H when the outer phase is 
collapsed and (j>H > 0. In this case the two coexisting 
phases are dense (Fig. ]^ . When = the outer phase 
is swollen and g{z) does not diverge at H (Fig. 0J)). A 
rough approximation yielding closed form expressions for 
g{z) for discontinuous brushes is described in Appendix 
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FIG. 8: H as a, function of a for different X2 values and TV 
3004^ 



THE COMPRESSION FORCE PROFILE OF 
A "TWO-STATE" BRUSH 



FIG. 9; The compression force profile for a brush with xi4') = 
1/2 + 1. 05(^^, a = 120 and iV = 300. The uncompressed brush 
is in a single phase state {cj>o < 4>- ) •''^ 



The surface force apparatus allows to measure the 
restoring force arising upon compression of a brush. For 
brushes of polymers characterized by a constant x the 
force increases smoothly with the compression and the 
force profile is essentially featureless. When the brush 
consists of polymers characterized by x(</') the comj 
sion can induce a vertical phase separation even if 
concentration profile of the brush is initially continu 
The existence of vertical phase separation, be it comj 
sion induced or not, gives rise to distinctive regime 
the force profile. In particular, the slope of the forc( 
distance curve in different compression regimes cai 
markedly different. In such experiments H is determ 
by the compressing surface rather than by a. Acc 
ingly, is set by the normalization condition (jSJ 
not by n((f>H) — 0. The compression increases Fchain 
the restoring force per area a is 



f{H) = 



dH 



In the following we obtain this force law for the cas 
compression by impenetrable, non-adsorbing surface 
For a brush with a continuous 4>{z) 



kT 



/oo(0) TTI' 



00 

2a2 8iV 



L 



z\<P)g{<P)^d<P+ 



b) 



C) as. 



chain 

kT 



obtained by invoking E{z',z)dz — z'^-k'^ /{8N). ] 
g{z) is given by while z(0) and dz/d(j> are spec 
by l(T3|) . f{H) is calculated numerically subject to 
constraint H16() . When the concentration at the wall 
exceeds 4>^ , the brush undergoes a vertical phase sep 
tion and 4'{z) is no longer continuous. In this case 
assumes the form 



8 10 12 14 



FIG. 10: (j}{z) plots corresponding to the three regimes in Fig. 
|5| (a) single phase swollen phase, H — 12.9, (b) a coexistence 
of a dense and a dilute phase, H = 6.1, (c) single dense phase, 
H = 3.7. In every case a = 120 and N = 300<^ 
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When the conditions permit a vertical phase separa- 
tion within the brush, it can take place in two ways. It 
can occur when the grafting density exceeds a certain 
critical value thus causing (j)Q > </>_. Alternatively, it 
can also take place as a result of compression when the 
grafting density does not lead to phase separation in the 
unperturbed brush. The development of (t>{z) and f{H) 
for this second case is depicted in figures El and ^| re- 
spectively. Initially, the brush retains the single phase 
structure and the associated force law. When the com- 
pression enforces (po > (j)^, a. vertical phase separation 
occurs and is signalled by a weaker slope of the f{H) vs. 
H curve. Stronger compression causes complete conver- 
sion to a dense phase thus causing an abrupt increase in 
fiH). 



VI. THE KARLSTROM MODEL AND THE 
DISTRIBUTION OF MONOMERIC STATES 

Thus far, our discussion concerned brushes character- 
ized by an arbitrary Xeffi4')- We now illustrate these 
considerations for the case of the Karlstrom modeli We 
focus on this model because of its simplicity and its semi- 
quantiative agreement with the phase diagram of aqueous 
solutions of PEO at atmospheric pressurai and the mea- 
sured x(2^j 4') ^ Within this model the monomers exist in 
two states: a polar, hydrophilic state (A) and an apolar, 
hydrophobic state (B). The two intercovert via internal 
rotations. In the N oo limit the corresponding Flory- 
type free energy is 



kT 



(l-0)ln(l-0) + 

(/)[plnp+ (1 -p) ln(l -p)] + 

(/)pAe 0(1 - (/)) [pxAS + (1 - p)XBs] + 

cf'xABPil-p) (31) 



where p denotes the fraction of monomers in the A state. 
The first term allows for the mixing entropy of the sol- 
vent while the second reflects the mixing entropy of the 
A and B states along the chain. Ae is the energy differ- 
ence between two states and 0pAe allows for the effect 
of interconversion between them. The fourth term is the 
generalization of the x4'{^ ~ 4') in ttio Flory free energy 
to allow for the interactions of the two states with the 
solvent S. The interactions between A and B states gives 
rise to the last term. 

The equilibrium value of p for a given is specified by 
dfoa/dp — leading to 



P 



1-p 



= cxp [-Ae - (1 - <j>){xAS - XBs)- 
cI)xab{1-2p)] 



(32) 



The parameters used by Karlstrom^ to fit phase diagram 
of PEO in water are xas = SO.O/T, xbs = 684.5/r, 
XAB = 155.6/T, Ae = -625.2/T+ln8. In the remainder 
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FIG. 11: A plot of the fraction of hydrophilic states, p, vs. z, 
in the Karlstrom model^ for different a. In every case xas = 
80.0/r, XBs = 684.5/r, XAB = 155.6/r, Ae = -625.2/r + 
InS, iV = 300 r = 60° Cp^ 



of this section we consider brushes at T = 60° C. The 
equilibrium Xeffi't') is obtained by equating H31(l and l|32(l 



XeffW 



PXAS + {i-p)XBS + 

-T—: [XABPil -p)+ pAe] 
1 - 

p Inp + (1 — p) ln(l — p) 



1 



where p(0) is specified by 132|l . x = ' 
with H32(l yield the equilibrium value 



(33) 



together 



X(0) = PXAS + (1 - P)XBS - XABP{1 - p) 



(34) 



Both Xeff{4>) and x(0) can be expanded in powers of 
(j). The coefficients in the expansion depend on the pa- 
rameters, Ae, XAS, Xbs, Xab- High order terms are of 
negligible importance and x((/i) ~ 0.48 4- O.Zlcf) + 0.07(/)^ 
provides a good approximation for x(0) • 

Equation H32I) allows to relate the volume fraction (f) to 
p, the fraction of hydrophilic A states, as 



In^ 



Ae 



XAS - XBS 



XAS - XBS - Xab{1 - 2p) 



(35) 



Accordingly, the exchange chemical potential can be 
specified in terms of p, i.e. A/i(0(p)). In turn, eq. I|13|) 
enables us to obtain z — z{p). To this end we invoke 
two boundary conditions: (i) In this range of parame- 
ters the brush is swollen and vanishes at the outer 
edge, 4>H = 4>{p = Ph) = where ph is the value of 
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FIG. 12: The overall ^(z) and the corresponding concentra- 
tion profiles of the monomeric states <j}A{z) — p{z)(f){z) and 
4>b{z) = [1 — p{z)](j){z) within the Karlstrom model for the 
conditions specified in Fig. 1111 — 

p at the height H . (ii) At the grafting surface we have 
A/x(z = 0,p = po) = BH^, where po is the value of 
p at z = 0. In addition we utihze the conservation of 
monomers as given by H16() . The corresponding plots of 
p = p{z) as well as the concentration profiles of the two 
states are depicted, for different u, in Fig. ^] and in 
Fig. El Since all brushes considered are swollen, with 
(/)_ff — 0, the p values at the outer edge of the brush, 
z ~ H are identical, p = pn- Increasing grafting den- 
sity leads to higher concentration at the grafting surface. 
This favors the hydrophobic B state and lower p at the 
surface. For the chosen parameters, the minimal value of 
p, corresponding to a PEO melt (0 = 1), is = 0.45. 

VII. DISCUSSION 

In this article we presented a common framework for 
the analysis of the structure of brushes of neutral water- 
soluble polymers (NWSP) that own their solubility to the 
formation of H-bonds with water molecules. Our anal- 
ysis concerned a family of two-state models developed 
for PEO but applicable, in principle, to other NWSP. 
The particular aspects of the models were grouped into 
Xef f {T, (p) thus allowing for a unified discussion of the 
brush structure within these models. Significant part 
of the discussion concerned brushes exhibiting vertical 
phase separation that can occur for polymers capable of 
a second type of phase separation. In particular, we ex- 
amined the distinctive behavior of plots of (z) and \/{z^ 
vs. a and the compression force profiles associated with 
such brushes, ^(z) and its moments are insensitive to 
the precise form of g{z) and the SCF analysis recovers 
the results obtained earlier^^ by the use of the Pincus 
approximation i i nSii^m In marked contrast, the compression 



force law does depend on g{z) and a full SCF analysis is 
necessary in order to obtain the correct results. These 
features are useful criteria for the occurrence of verti- 
cal phase separation. Such criteria are of interest be- 
cause of indirect experimental indications that brushes 
of PNIPAM exhibit this effect. Among these indica- 
tions the following are especially noteworthy. First, is 
an early study by Zhu and Napper^^ of the collapse of 
PNIPAM brushes grafted to latex particles immersed in 
water. This revealed a collapse involving two stages. 
An "early collapse", took place below 30° C, at "bet- 
ter than ^-conditions" , and did not result in flocculation 
of the neutral particles. Upon raising the temperature 
to "worse than 6'-conditions" the collapse induced floc- 
culation. This indicates that the colloidal stabilization 
imparted by the PNIPAM brushes survives the early col- 
lapse. It lead to the interpretation of the effect in terms 
of a vertical phase separation within the brush associated 
with a second type of phase separation as predicted by 
the n-cluster model. More recently, Balamurugan et al^ 
used Surface Plasmon Resonance (SPR) and contact an- 
gle measurements with water to characterize PNIPAM 
brushes "grafted from" a self-assembled monolayer on 
gold. The brush properties were studied at between 10° 
C and 40° C. The SPR measurements indicated gradual 
variation of the brush properties. In marked contrast, the 
contact angle measurements revealed an abrupt change 
at ^ 32° C. Both experiments are consistent with a brush 
undergoing vertical phase separation such that the outer 
phase is hydrophilic. Within this picture, the abrupt 
changes reported in the wetting and aggregation behav- 
ior correspond than to the onset of a single, dense and 
hydrophobic phase. This picture is also supported by re- 
cent study of the phase behavior of PNIPAM by Afroze 
et al& Early study of the phase behavior of PNIPAM 
in water, by Heskins and Guilletf2M, identified a LOST at 
(j)c — 0.16 and Tc — 31.0° C. In marked contrast, the 
work of Afroze et al^'^ identified PNIPAM as a polymer 
undergoing a second type of phase separation. In partic- 
ular: (i) While the LOST of PNIPAM depends on N the 
LOST occurs around ^ 27 - 28° C and <j)c ^ 0.43 (ii) 
In the limit of 0, the phase separation occurs, de- 
pending on N, between 30° C and 34° C. Thus, the phase 
diagram of Afroze et al suggests that a vertical phase sep- 
aration is indeed expected in brushes of PNIPAM. Sys- 
tematic Neutron Reflectometry (NR) studies of PNIPAM 
brushes will eventually provide clearer picture of the sit- 
uation. Early studies were hampered by high polydisper- 
sity as well as difficulties in determining N and a.^^ With 
this in mind, the NR results revealed that the structure 
of PNIPAM brushes in acetone was very different from 
their structure in water, both at 20° C and at 55° C. 
In acetone the concentration profile was smoothly decay- 
ing while in water it consisted of a narrow, inner, dense 
region and an outer, extended and dilute region. More 
recent work utilized NR,'*" to study samples with lower 
polydispersity and higher grafting density between 20° C 
and 40° C. Importantly, the results indicate that most of 



10 



the conformational change occurred between 28° C and 
34° C though the corresponding concentration profiles 
were not reported. The results suggest a repartitioning 
of the monomers between a dilute outer tail and an inner 
dense region.^- Finally, recent experimental results of Hu 
at al^"^ are suggestive of the predictions obtained above 
concerning the variation of (z) upon decreasing a (Fig. 
El). Hu al studied the thickness of a PNIPAM brush 
grafted to spherical microgels of copolymers of PNIPAM 
and acrylic acid 2-hydroxyethyl ester (HEA). The mi- 
crogels shrink as T increases from 24° C to 36° C thus 
inducing a decrease in a. Remarkably, the thickness of 
the brush initially decreases in the range 27° C to 32° C 
but subsequently increases upon further heating in the 
range 27°-32° C. Unfortunately, in this experiment it is 
impossible to separate the effects due to change in T from 
those to the change in a because the T is used to tune a. 

In confronting experimental results with theoretical 
predictions concerning the vertical phase separation it is 
important to note two points. First, polydispersity in N 
and in a can smooth out the discontinuity in 0(2:). These 
factors give rise to domains with different 0o. In turn, 
the altitude of the discontinuity in these domains will 
also differ. NR measurements will reflect the weighted 
average of these domains. Polydispersity in N can lead 
to further strengthening of this effect if Xef / depends on 
N . Second, the vertical phase separation involves a first 
order phase transition. Accordingly, nucleation dynam- 
ics may play a role and it is necessary to allow for the 
possibility of long relaxation times. 



2B 



<t>+(Ht) 



A^{(l))d(l) 



2B 



The evaluation of (|22|l in the case of a continuous (/)(z) 
involves introducing the variable A/i = B{H^ ~ z^) ^nd 
invoking z = \J H"^ — A/i/_B. Integration by parts leads 
to 



cr 1 

iV^2B 
a 1 



z(j)[z)dlS.pL 



Afi=BH-' 



A^(i(z0(z)) 



lAfi=0 

a 1 /■'^V , 5z 



Na^ 2B 



ct>H 



d(j) 



z ) A/ifi0 (A3) 



Using 



dz 



d(p B^m - A^i/B 
we obtain the final expression 



(A4) 



I r<t>a 0A/i 



9Ap 



Na^ 2B2 7o ^i?2 - A^i/S 
1 



2B 



Afiy/H^ - Afj./Bd4> (A5) 



APPENDIX A: CALCULATION OF THE 
AVERAGE THICKNESS 

For a (j>{z) 1211) is evaluated using integration by parts 
utilizing dAjj, = —2Bzdz: 



a 1 



Na^ 2B Jo 



EH-' 



{z)dA^jL 



For a discontinuous (\){z) (|21ll this procedure leads to 



z(j){z)dz + / z(f){z)dz 



a 1 

(\)(z)dA^i 



(z)dA/^- 



^[H --{H -H,) 



where A/z is specified by (fnj) . 

Following a similar procedure for a discontinuous 4>{z) 
(EH) yields 



G 1 



Ht r-H 

z'^<l)[z)dz -\- I z'^(f){z)dz 

JHt 



BH^ 



Na^ 2B {JsiH^-H! 

z(f){z)dAii 



zif){z)dAfi~ 



(t)Afi 



2i?2 - Afi/B 



— I Afi^m - A^l/Bdc^ 



<l'+{Ht) 

1 r^-iH.) 0Am^ 



d<j)^ 



2^2 7o ^/W^AjI/B 



Afiy/H^ - A^i/Bd(t) 
-Ht {H^ - Hf) 



(A6) 



11 



APPENDIX B: CALCULATION OF THE 
DISTRIBUTION OF ENDS FUNCTION 

Upon introducing the variables p — — z^, t = — 
z'^ and g{z')dz' = —f{t)dt eq. Q assumes the form of 
an Abel integral equation 



v{p') 



and its value at p' = is 



f-it)dt 



(B7) 



v{p) = 



f{t)dt 



(Bl) 



where v{p) = 



2Nw 



) whose solution IBlll is 



V VP Jo VP-t dt 
It is convenient to rewrite the integral in HB2() as 



(B2) 



1 



dv{t) 



Vp- i d,t 



dt = 



viz) 



dv 



v{H) y^p-t{v) 



2Na^ J^^ ^A/i(0) - A^(0') 



(B3) 



where Ap{(j)) is given by H13|l . Substituting this integral 
into (|B2|I while noting that v{0) — ^ Jj'a i 4>h and g{z) = 
2zf{p), leads to g{z) in the form 123|l . This equation 
yields a simple form z{(f)) rather than for (j^iz). Thus, it 
is naturally to consider g{z) as a parametric function 



_ AM(0)/i? 



z{(t>) 



B 



(B4) 
(B5) 



These two equations are augmented by (|16|l relating a to 

00- 

In the case of discontinuous (j){z) it is necessary to ob- 
tain g{z) in the outer and inner phases separately. For the 
outer phase the introduction of the variables p = H'^ — z'^, 
t = H'^ - z'2 and g^{z')dz' ^ ~f^[t)dt transforms ^ 
into an Abel integral equation IjBip whose solution is 

To obtain g+{z) at the inner phase we substitute H26|l 
into H25|) and transform the first term of H25|) into an Abel 
integral ljBl|) by introducing the variables p' = Hf — z^, 
t' = - z'2 and g+{z')dz' = -f+{t)dt. This leads to 



TTcr 
2iVa3 

TTCT 
2N^ 



-H^-t 



(B8) 



where eq. (|B1|I with p = — Hf was used in order to 
calculate the second term in (IB8II. In addition 



dvjp') 
dp' 



TTCT C?(/)+(z) 

2N^~df/~ 



f-it)dt 



{H^ - +t' -t) 



3/2 



(B9) 



All together, upon substituting l|B7|) . (|B8|) and l|B9|) 
into (IB2II we obtain 



27Va3 



1 



dt' 

Vp^Jo 



27r 



!-{t)dt 



[m-H^ + t'-t) 



In the first integral we express p' and in terms of 
(|13|l and we change the order of integration in the second 
integral obtaining 



Uip) = 



1 

2^ JO 
p' 



f-(t)dtX 
dt' 



lo (iJ2 _ ij2 ^ ^, _ ^)3/2 

The inner integral in (|B11|1 is 



(Bll) 



f-{t)dt _ fP f+{t')dt' 



Vp^ 



f, (B6) 



where f-{t') = g_(z')/2z', and p ^ p' 
v{p') in the solution of the Abel equation (|B1|I is 



H^-Hf. Thus, 



Vp' - t' {H^ - Hf +t' -t) 
2V7 



3/2 



(B12) 
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leading to 



f+ip') = 



2Na^ 



.{Ht) - 4>-{Ht) 



+ - 



f^{t)dt 



(p' + m- t) H'f- 1 



Substitution of f-{t) = g_(z)/2z from H26|l yields the 
final expression for g+{z) in the form 



9+iz) 



-{Ht 



TT 



4,'Jt) 



Hf - z2 X 



0« VA/i(^0 - AAi((/.:.') 



(B14) 



It is convenient to express g+{z) in terms of the con- 
centration of the dense phase, 0+ utilizing z = 
H'^ — A/i((/)+)/_B. Introducing the variables t = 
Afi{(j)'J/B and = 1/ B{dAfi{(j)'J/d<j)'Jd(l}'_ in the last 
integral leads to (|27|l . 



APPENDIX C: MODEL CONCENTRATION 
PROFILES 

A rough approximation allows to obtain analytical ex- 
pression for g{z) of a brush in the presence of a vertical 
phase separation. In the dense phase the variation of (j) 
is slow and we can approximate it as constant, = 4>o- 
In the outer phase dilute phase is rather low and we 
can neglect nonconstant terms in x{(j>)- Thus 



(j>o,0< z <Ht 
{z),Ht<z<H 



(CI) 



where 4){z) is determined by 1)1 8|l with x — const, while 
the value of 0o is set by©- As before, g{z) in the outer 
phase is determined by 12()|) while in the inner phase it is 
determined by ljB14|l . 

First, consider the case of x = 1/2 leading to Ap((/)) = 
- ln(l - 0) - « (jP/2 and eq. ^ yields 



'2B 



(C2) 
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FIG. 13: Comparison between the exact (j>[z) and g{z) and 
their approximate values as calculated from for the case of 
X(0) = 1/2 + 1.05<^^ a = 120 and N = 3004^ 



Substitution of (|C2|I into ljB14|l gives the expression for 
g{z) in the inner phase 



za 



'2B arctan - 



(C3) 



where <p-{Ht) = y/2B{H^ - Hf) and (t)+{Ht) = <j)o. 
This function is discontinuous and diverges at the phase 
boundary z ^ Hf . The value of 0o can be found from Q 



arcsm — 

Ht \2 H 



m Hf 



(C4) 

When X = eq. H26|) specifies g{z) at the outer phase. 



Ht<z<H 



(C5) 



and eq. I|B14|) for g{z) in the inner phase, < z < Ht, 
yields 



7Va3 



b+{Ht) - ^-(Ht) 
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Again, g(z) of the inner phase diverges at the phase 
boundary. The performance of this approximation is il- 



lustrated in Fig^] It captures the main features of 0(z) 
and the behavior of g{z) in the inner region. However 
g{z) at the outer region increases rather than decrease. 
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